Preparing a mechanical oscillator in non-Gaussian quantum states 
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We propose a protocol for coherently transferring non-Gaussian quantum states from optical field 
to a mechanical oscillator. The open quantum dynamics and continuous-measurement process, 
which can not be treated by the stochastic-master-equation formalism, are studied by a new path- 
integral-based approach. We obtain an elegant relation between the quantum state of the mechanical 
oscillator and that of the optical field, which is valid for general linear quantum dynamics. We 
demonstrate the experimental feasibility of such protocol by considering the cases of both large- 
scale gravitational-wave detectors and small-scale cavity-assisted optomechanical devices. 
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Introduction. — It is becoming experimentally possible 
to prepare a macroscopic mechanical oscillator near its 
quantum ground state by either active feedback or pas- 
sive cooling in optomechanical devices [ij. This activity 
has been motivated by (i) the necessity to increase the 
sensitivity of high-precision measurements with mechan- 
ical test bodies up to and beyond the Standard Quan- 
tum Limit (SQL) [2|, and (ii) the test and interpretation 
of quantum theory, when macroscopic degrees of free- 
dom are involved. However, for unequivocal evidences 
of quantum behavior, merely achieving quantum ground 
state, or preparing coherent /squeezed states, or overcom- 
ing the SQL is insufficient. In these situations, the os- 
cillator initially occupies a Gaussian state and remains 
Gaussian, and therefore its Wigner function is positive 
and can always be interpreted in terms of a classical 
probability. A true demonstration of the quantum be- 
havior requires non- Gaussian quantum states or nonlin- 
ear measurements [3, 0] • natural approach is to create 
nonlinear coupling between a mechanical oscillator and 
external degrees of freedom, e.g., probing mechanical en- 
I, coupling the oscillator to a qubit 
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(low) cavity photon number |12h14| . For optomechanical 

devices, this generally requires zero-point uncertainty of 

the oscillator displacement Xq to be comparable to the 

cavity linear dynamical range which is characterized by 

the optical wavelength A divided by the finesse J-', i.e.. 
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Since A ^ lO^^m and T < 10^, wc have Xq > 10~"'m, 
which is several orders of magnitude above the current 
technology ability. 

In this article, we propose a protocol for preparations 
of non-Gaussian quantum states which does not require 
nonlinear optomechanical coupling. The idea is to inject 
a non-Gaussian optical state, e.g., a single-photon pulse 
15l4l7j. into the optomechanical 



created by cavity QED 
devices. Possible configurations are shown schematically 
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FIG. f: (Color online) Possible schemes for preparing non- 
Gaussian quantum states of mechanical oscillators. The 
left is a Michelson interferometer, similar to an advanced 
gravitational-wave detector with kg-scale suspended test 
masses [l^ . figl ] . The right panel shows a small coupled-cavity 
scheme with a ng-scale membrane inside a high-finesse cavity 
01 • In both cases, a non-Gaussian optical state (a photon 
pulse) is injected into the dark port of the interferometer (lo- 
cal oscillator light for homodyne detection is not shown). 



in Fig. [TJ The radiation pressure induced by the photon 
pulse is coherently amplified by the classical pumping 
at the bright port, and the qualitative requirement for 
preparing a non-Gaussian state is 



X/iTXq) 
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Here = IoT/(huo) (/o the pumping laser power 
and LdQ the frequency) is the number of pumping pho- 
tons within the duration r of the single-photon pulse. 



and we gain a significant factor of 



compared with 



Eq. ([T]), which makes it experimentally achievable. This 
radiation-pressure-mediated optomechanical coup ling is 
similar to what was considered in Refs. How- 
ever, there are significant differences: (i) This protocol 
includes both finite interaction time and photon shape, in 
which case neither the rotatin g-w ave approximation [20[ 
nor the three-mode approach [2l| applies; (ii) To better 
model an actual experiment, we consider a continuous 
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measurement process rather than a single measurement 
at some given instant as assumed in Ref. . This takes 
into account all the information of the oscillator motion 
that is distributed in the output field, and thus allows us 
to prepare a nearly pure non-Gaussian quantum state of 
the oscillator; (iii) There are non-trivial quantum corre- 
lations at different times (non-Markovianity) due to the 
finite-duration photon pulse, which cannot be treated by 
the conventional stochastic-master- equation (SME) ap- 
proach |23l - l27l |. Here we develop a path- integral-based 
approach, and it applies to general linear quantum dy- 
namics and continuous measurement process. 

A simple case. — To illustrate the non-Gaussian state- 
preparation procedure, we first make an order-of- 
magnitude estimate of experimental requirements by con- 
sidering a simple case where the cavity decay is much 
faster than all other time scales and the oscillator can be 
approximated as a free mass. The corresponding input- 
output relations, in the Heisenberg picture, simply read: 

i{t)^m/m, Mt) = c^a,{t)+Fth{t), (3) 
6i(t)=ai(t), b2(t)=d2{t) + {a/h)x{t). (4) 

Here x and p are position and momentum; the coupling 
constant a = Sy/2{T / X) ^/hlojujo] ai_2 and 61^2 are in- 
put and output optical amplitude and phase quadratures, 
with [ai{t),a2{t')] = [bi{t),b2{t')] =iS{t-t'); aai is the 
back-action noise; Fth is the force thermal noise. 

Suppose at t = — r the oscillator was prepared in some 
initial Gaussian state \ipm) = J^^'4'm{x)\x) dx (the pro- 
cedure is detailed in Ref. (l^l). Subsequently, a photon 
pulse is injected into the dark port of the interferometer 
and starts to interact with the oscillator. During this 
interaction, phase quadrature 62(0 is continuously mea- 
sured by a homodyne detection, until the photon pulse 
ends at t = 0. If photon pulse (i.e., t) is short such that 
oscillator position almost does not change, we obtain: 

1(0) = X{-t), P(0) = P(-r) + K ii + Ah , (5) 

B2=A2 + kX{Q). (6) 

We have normalized the oscillator position and mo- 
mentum by their zero-point uncertainties: X = x/xq 
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{pq = ^Jhn-iuJml'2\\ 
which has an unccr- 



^J%dtd,{t){j = 1,2) 
tainty of unity (i.e., AAj=l): Bj = \/1/t J^^ dt bj (t) ; 

Ah = /"^ dtFth{t)/pq] K = ay^/h = 8V2y^TXq/X. 

Eqs. ^ and ([6]) describe the joint evolution of the os- 
cillator, the optical field and heat bath in the Heisen- 
berg picture (with Bj viewed as the evolved versions 
of Aj). They transform back into an evolution oper- 
ator of [/ = exp[i{KAiX -\- Ah-^)] in the Schrodinger 
picture. The corresponding density matrix of the sys- 
tem at t = is given by p = t/|'!/'o)|V'm)/5th(V'm|(V'o|C^'l', 
where \ipo) = /^^/'o(^2)|^2) c^^2 is the initial non- 
Gaussian optical state, and pth describes the heat 
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FIG. 2: (Color online) A schematic of the non-Gaussian state- 
preparation process. The interaction entangles the oscillator 
state and the optical state (depicted by their Wigner func- 
tions). Subsequent measurements of the optical fields dis- 
entangle the system and projects the oscillator into a non- 
Gaussian conditional quantum state. 



bath associated with Pth- Given homodyne detection 
of B2 with a precise result y, the oscillator is pro- 
jected into the following conditional state: Pm{y) = 



Tr 



th 



In the ideal case 



{y\U\iJo)\^m)Mi^m\mU^y) 

of negligible thermal noise, the conditional wave function 
'0^j(a;) of the mechanical oscillator is simply 



V'm(a;) = 4'oi.V - Kx)^m{x) 



(7) 



— the optical state is mapped onto the mechanical oscil- 
lator as illustrated in Fig. [51 A complete mapping oc- 
curs when ipm{x) « const, and this requires the momen- 
tum fluctuation due to optomechanical coupling be larger 
than the initial one, namely, k > 1 or equivalently 
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which justifies Eq. 

When thermal noise is considered, non-Gaussianity can 
still remain, as long as thermal noise induces a smaller 
momentum fluctuation than the optomechanical interac- 
tion. This condition, in the high-temperature limit — 
(Ah(t)Ah(i')> = ^m-f^kBT5{t - t')., reads 



A/ {TXq)yJ nth/ Qmy/i^mT < 8V2y^ 



(9) 



with Qm = '^m/lm thc mcchauical quality factor and 
rith = ksT / {huim) the thermal occupation number. 
These two conditions set the benchmarks for a successful 
non- Gaussian state-preparation experiment. They can 
be satisfled with experimentally feasible specifications as 
shown in Table HI in which the first row is similar to the 
case of large-scale gravitational-wave detectors [l^ and 

TABLE I: Possible experimental specifications 
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large scale 


l^m 


6000 


4kg 


IHz 


10* 300K 1ms 


small scale 


Ifxrci 


10* 


Ing 


lO^Hz 


10^ 4K 0.01ms 
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the second row is for small-scale optomechanical devices 
(e.g., the one in Ref. 0). These qualitative results will 
be justified by a rigorous treatment below. 

General formalism. — In general, the optomechanical 
interaction strength is finite and the oscillator has non- 
negligible displacement during the interaction, the cav- 
ity bandwidth can be comparable to the mechanical fre- 
quency, and thermal noises can be non-Markovian. All 
these factors obstruct finding a finite set of variables sim- 
ilar to {X, P, Ai,A2) that satisfy a closed set of equations 
[cf. Eqs. ([5]) and ([H])]. It is therefore hard to determine, 
a priori, the finite number of observables that one has 
to measure to project the oscillator into a desired condi- 
tional state. 

To address these issues, we adopt the Heisenberg pic- 
ture starting from t = — oo, and write down the initial 
density matrix as pin = /5,„(— oo) (8) Po (8i Pth- Details of 
Pm{—oo) for the oscillator and whether the initial state is 
truly a direct product, do not matter, because the system 
is stable, and the initial position and momentum will de- 
cay away after several mechanical relaxation times. For 
the optical state, we consider an arbitrary spatial mode 
given by /(x/c), whose annihilation operator is 



The generating function for the oscillator state is then 



t^f_^dtfit)[a,{t)+za,it)]/V2. 



(10) 



A general state of this mode can be written in the P- 
representation a.s po — J c?C^(C)IC)(C|j where vector ^ = 
(3?[C],5[C]) and |C) = cxp[Cft - Cf]|0). 

A continuous measurement of the output optical 
quadrature y{t) = cos9bi{t) + sinO b2{t) for t £ (— oo,0], 
projects the entire system into a conditional state: 



Pc[y{t)] = Vy PinVy/Tr[Vy PinVy]. 



(11) 



The operator Vy projects the output field into the sub- 
space where y{t) agrees exactly with the measured re- 
sults y{i). To simplify output correlations at different 
times, we can causally whiten y{t) into z{t) such that 
{z(t)z(t')) = S{t - t'), as detailed in Ref. i28|. Since the 
output quadratures at different times also commute, i.e., 
[z{t), z{t')] = 0, the projection Vy can then be expressed 
as the product of Dirac-(5 functions that project each z(t) 
into its measured value z{t): 



cx)<t<0 

v[i]c^^[if^^dtm[m - -At)]} ■ (12) 

with j 2?[^] denoting the path integral. This allows us to 
take the entire measurement history for z (or equivalently 
y) and project into the corresponding subspace in a single 
step, instead of having to successively project output- 
field degrees of freedom continuously at each time step 
as in the case of SME approach, thereby allowing a non- 
Mar kvonian input field. 



^[a;z(i)] =Tr 



e'"""Pc[^(i)] 



(13) 



where a. = {ax,ap), Xq = {x(0),p{0)), and superscript ' 
denotes transpose. From Eqs. pT|) and p^ . we have 

j dCP{C) Jv[^]e'^'^'^~'^^'^'^^0\e'^\0} (14) 

with B = ax'„ + J°^dt^{t)[z{t) - z{t)]. This can be 

evaluated by decomposing Xq = R + j-oo'^^-^^^^)^^^) 
where K = [Kx, Kp) are causal Wiener filters, K{—t) = 
(O|z(i)xo|0) and R = {Rx,Rp) are parts of displace- 
ment and momentum uncorrelated with the output: 
{0\Rx,pz\0) = 0. Completing path integral, we obtain 

J=J d(^e-l°''f-='^' + \\''-'^CL'f]/2+ta{CY+C-(^+x'J p(^)(i5) 

HereljajpEE J_^a(t)a*{t)dt and we have defined vectors 

7 EE [r,R] and L = (5R[L],3[L]) with L{t) = [f,z(t)], 
which characterize the extent of photon mode influence 
on the fluctuations of x{0) and p(0), and output field 
z; quantities Vc = {0\BFR\0) and Xc = {xc,Pc) = 
j'^^dtK{—t)z{t) arc the conditional covariance matrix 
and means of x{0) and p{0) when the optical state is vac- 
uum. The resulting conditional Wigner function reads 



W[x;z(t)] ^ / dCe 



,-[xV7'x' + \\z-2CL'f]/2 



PiC) (16) 



with X = a; — — C*7 ^ C7*- This formula directly re- 
lates the injected optical state to the state of the mechan- 
ical oscillator. In deriving it, we only use the linearity of 
quantum dynamics rather than specific equations of mo- 
tion. For cavity-assisted optomechanical system, one can 
obtain -7, Vc, K and L from input-output relations in 
Refs. [29I43H by using formalism developed in Ref. [2^. 

Single-photon case. — As an example, we consider the 
simplest case of a single-photon injection, with po = 
|1)(1| and P(C) = eJ^^^d^S^^^O/dCdC- From Eqs. (HH), 
it gives 



W[x;z{t)] = 



l-jY-'ji -\\Lf + \jV-^Sx' + Z\ 



1 



27rVder 



: exp 



--6xY-^Sx' 
2 



(17) 



where 6x = x — Xc and Z = j'^^dt z(t)L{t). This 
Wigner function depends on the measurement result z{t), 
t G (—00, 0] through four quantities, the two components 
of Xc (through 5x) and the real and imaginary parts of 
Z: Z determines the shape of W , and x^ describes the 
translation of W. The random vector Z = (5R[Z],S[Z]) 
has a two-dimensional probability density of 



w[Z] 



1 



\L\ 



ZZ' 



27rVdet Vl 



cxp[-ZV^iZ72], (18) 
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reconstruction, which is crucial for reveahng those nega- 




tivity regions. 
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FIG. 3: (Color online) Distributions of measurement results 
(left panels) and the corresponding Wigner function of the os- 
cillator given the most probable measurement results (middle 
panels) and less probable results but with a significant non- 
Gaussianity (right panels). The upper panels show the case 
of non-Gaussian state-preparation with future gravitational- 
wave detectors, and the lower panels for small-scale devices. 
We used normalized coordinates (with respect to Xq and pq) 
and introduced Qg = ^hm/a^. 



where matrix Vr = f'^ dtL'L. 

The pre-factor in the Wigner function [cf. Eq. pT)) ] is 
a second-order polynomial in a;, which resembles that of 
a single-photon. For strong non-Gaussianity, significant 
7 and (making 7 terms in the pre-factor to prevail) 
are essential — these physically correspond to requiring 
that the photon mode must influence the fluctuation of 
X and p, as well as z strongly. It in turn requires the 
photon coherence time to be comparable to the measure- 
ment time scale characterized hy K. It is possible for 
small-scale optomechanical devices with high-frequency 
mechanical oscillators. The correspondi ng p hoton can 
be generated by a cavity QED scheme IB" 12 1- While for 
large-scale gravitational-wave detectors, the time scale is 
~ 10 ms and it is challenging to create photons with 
comparable coherent length. However, developments of 
low-frequency squeezing source (s^ will eventually solve 
this issue. 

With Eq. (fT7|) . we can justify the simple-case qualita- 
tive results. We use the same specifications listed in Ta- 
blelH As an example, we assume a photon shape of /(i) = 
y2^g(7jH-ja;j)t g^j^^ specify that w//27r = 7//27r 70 Hz 
in the case of future gravitational-wave detectors, and 
LOf/uJm = 0.1, ^f/uim = 0.3 for small-scale experiments. 
The Wigner functions for some given measurement re- 
sults are shown in Fig. |3l In both cases, there are neg- 
ative regions in the Wigner function, which is a unique 
feature of the quantumness. The prepared non-Gaussian 
quantum state can be independently verified using the 
quantum tomography protocol developed in Ref. [s^ 
that allows sub-Heisenberg accuracy of Wigner function 
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